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Abstract: In this paper, we show that approximate nonlinear self-adjointness of per- 
turbed PDEs is equivalent to nonlinear self-adjointness of the corresponding unperturbed 
PDEs, which provides a shortcut to discriminate approximate nonlinear self-adjointness. 
Moreover, a succinct approximate conservation law formula by virtue of the known con- 
servation law of the unperturbed PDEs is given in a explicit form. We consider several 
examples to illustrate the results and to show that approximate nonlinear self-adjointness 
can generate new approximate conservation laws. 
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fri 1 Introduction 

The discovery of conservation laws played a pivotal role in the comprehensive study of the 
^ ■ properties and solutions of partial differential equations (PDEs) arising in many disciplines 

of the applied sciences including physical chemistry, fluid mechanics and etc. For example, 
the knowledge of conservation laws is useful in the numerical integration of PDEs such as 
^ ■ controlling numerical errors. An infinite sequence of conserved densities is a predictor of 

the existence of solitons and complete integrability which means that the PDEs can be 
solved with the inverse scattering transform method pp. Hence, in order to achieve these 
goals, one needs efficient methods to compute conservation laws. 

The relationship between symmetry and conservation law provides a feasible method 
to construct conservation law of PDEs. The celebrated Noether theorem gave a systematic 
way to determine conservation laws for system of Euler-Lagrange equations which relied 
on the availability of a Lagrangian [2]. Later, new methods were proposed which pose 
less restrictions on the knowledge of Lagrangian. Kara and Mahomed presented a par- 
tial Noether approach for Euler-Lagrange type equations [5J. In [2111], Anco and Bluman 
proposed multiplier method to construct local conservation laws for PDEs expressed in a 
standard Cauchy-Kovalevskaya form, which did not count on the existence of Lagrangian. 
Recently, Ibragimov suggested the general concept of nonlinear self-adjointness for con- 
structing conservation laws which embraces strict self-adjointness, qusi self-adjointness 
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and weak self-adjointness stated earlier and these methods have exerted effectiveness 
when applying to some physical PDEs [6Tfl2]. 

On the other hand, combination of symmetry group theory and perturbation analysis 
techniques leads to emergence of approximate symmetry which aims to study perturbed 
PDEs possessing few exact symmetries [THUS], thus it is natural to extend the above 
procedures for usual PDEs to handle approximate conservation laws associated with ap- 
proximate symmetry. In [16], Kara, Mahomed and Unal studied how to construct approx- 
imate conservation laws for perturbed PDEs via approximate Lie-Backlund symmetries. 
In [T7j, a basis of approximate conservation laws for perturbed PDEs was discussed. John- 
pillai et.al [18] showed how to construct approximate conservation laws of approximate 
Euler-type equations via approximate Noether type symmetry operators associated with 
partial Lagrangian. Quite recently, we introduced the notion of approximate nonlinear 
self-adjointness and studied its properties to construct approximate conservation laws 
which extends the results in unperturbed case [T9] . 

The purpose of this paper is to study how to discriminate approximate nonlinear 
self-adjointness of perturbed PDEs besides the definition. We show that approximate 
nonlinear self-adjointness of perturbed PDEs is equivalent to nonlinear self-adjointness of 
the corresponding unperturbed PDEs. We also give a specific approximate conservation 
law formula by virtue of the known conservation law of the unperturbed PDEs. 

The outline of the paper is arranged as follows. In Section 2, some related basic 
notions and principles are reviewed and the main results are given. In Section 3, several 
examples are considered to illustrate the results. The last section contains a conclusion 
and discussion. 

2 Main results 

In this section, we first recall some notions and principles related with approximate non- 
linear self-adjointness, and then give the main results of the paper. 

2.1 Preliminaries 
2.1.1 Notations 

We first define some notations which will be used throughout the paper. 

Let x = (x 1 , . . . , x n ) be independent variables and u = (u 1 , . . . , u m ), v = (v 1 , . . . , v m ) 
and w = (w 1 , . . . , w m ) as three dependent variable sets, u) = . . . , m ), if>(x, u) = 
(ip 1 , . . . , ijj m ) and (fi(x, u) = ((p 1 , . . . , (f m ) be three m-dimensional vector functions, whose 
different order derivatives are denoted as follows 

A (1) = {Af 1 }, A (2) = {Af lfa }, A (r) = 

where A£_ 4 = D^D^ . . . D ia (A a ') and a = 1,2, ... ,m. The symbol "A" denotes the 
dependent variables u,v,w and vector functions (f),ip,(p. Hereinafter, Di denotes the 
total derivative operators with respect to x l . For example, for one dependent variable 



u = u(x, t) with x = x 1 , t = x 2 , one has 

d d d d 

D x = — — h u x — + u xx - h u xt - h 

ox ou ou x ou t 

2.1.2 Approximate nonlinear self-adjointness 

Consider a system of m PDEs with rth-order 

E a = El(x, u, it(i), • • • , tt (r) ) + eJ5„(x, u, it (1) , • • • , w (r) ) = 0, (1) 

where e is a small parameter and a = 1,2, ... ,m. System ([1]) is called perturbed PDEs 
while the system which do not contain the perturbed term eE\(x, u, U(i), • ■ ■ , U( r )) in the 
form 

El{x,U,U(x), ■ ■ ■ ,U(r)) = 0, (2) 

is named by unperturbed PDEs. Note that the summation convention for repeated indices 
is used throughout the paper if no special notations are added. 
Let L be the formal Lagrangian of system (T5]) given by 

C = vPEp(x,U,U {1 ),--- ,U(r)), (3) 

then the adjoint equations of system §2§ are defined by 

{El)*{x,u,v,u {1) ,v {1) , - ■ ■ ,u (r) ,v {r) ) = — = 0, (4) 
where 5/5u a is Euler-Lagrange operator written as 

d - + Z ( -iy Dh ...n°. (5) 



5u a Q u a ^ > 'i ^Q u a 
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Definition 2.1 (Nonlinear self-adjointness JEj) The system (TJ|) is said to be nonlinearly 
self- adjoint if the adjoint system Op is satisfied for all solutions u of system (TJ|) upon a 
substitution 

v a = <f a (x, u), a — l,...,m, (6) 

such that <f(x, u) ^ 0. 

Here, ip(x,u) ^ means that not all elements of (f(x,u) equal zero. Definition 12.11 
is also equivalent to the following identities holding for the undetermined coefficients (or 
functions) 

(El)*(x,U,V,U m ,V (1) ,--- ,U(r),V( r ))\ {v= v,..., V[r)=nr)} = *i*%(x, U, U {1) , ■ ■ ■ ,U (r) ), (7) 



or equivalent form [5{yP E®) / 5u a }\ {v=v v =v } = A^i?°, which is applicable in the com- 
putations. 

In [19], we extended the above procedures to deal with perturbed PDEs. The formal 
Lagrangian C of perturbed system (CD) is given by 

Z = w p [E%(x, u,up), ■ ■ ■ ,u {r) ) + eE^(x,u,u {1) , ■ ■ ■ ,u (r) )], 

then the adjoint equations of system ([1]) are written as 

E* a {x, u, w, it(i), w {1) , . . . , w (r) , w {r) ) = — = 0. (8) 

Definition 2.2 (Approximate nonlinear self-adjointness JTty) The perturbed system (QP 
is called approximate nonlinear self-adjointness if the adjoint system (G|) is approximate 
satisfied for all solutions u of system (QJ) upon a substitution 

w a = ip a (x, u) + e(p a (x,u), a=l,...,m, (9) 

such that not all i\) a and 0°" are identically equal to zero. 

Definition l2.2l extends the results of unperturbed case and is equivalent to the following 
identity with undetermined parameters (or functions) A^ and fj,^, 

E* a (x, u, w, u {1) , w (1) , . . . , u (r) , w {r) ) l{w=i>+e . ^ w ^ (r)+£ , (r)} 

- [(A£ + e(iZ)E%(x, u, u {1) , . . . , u (r) ) + eAjSj(x, u, u {1) , . . . , u (p) )] = 0(e 2 ), (10) 

which provides a computable equality for approximate nonlinear self-adjointness of per- 
turbed PDEs. 

Definition 2.3 (Approximate conservation law WU$) The vector T = (T 1 ,T 2 , . . . T n ) de- 
fined by T l = T ' + eT{ + ■ ■ ■ + e k T' k is an approximate conserved vector of system (QJ) if 
T % satisfies approximate equation D{T l = 0(e k+1 ) for all solutions of system (T7]) ; which 
defines a kth-order approximate conservation law. 

The major consideration of the paper concentrates on the first-order approximate 
conservation law which is defined as D{T l = 0(e 2 ) with T l = Tq + eT[. 

The following theorem will be used to construct conservation laws for both unper- 
turbed and perturbed cases [2Tj . 

Theorem 2.4 Any infinitesimal symmetry (Lie point, Lie-Bdcklund, nonlocal) 

d , d 



X = C(x, u, u {1) , . . . )— + rf{x, u, u {1) , ■■■) duC 



of system (TJ|) leads to a conservation law Di{C l ) = constructed by the formula 
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du° 



+ ..., (n; 



where W a = rf — ^u c „- and K, is the formal Lagrangian. 



Generally speaking, the term £*/C with /C in the form ([3]) can be omitted because it 
vanishes identically on the solution manifold of the studying PDEs. 

The following properties about approximate nonlinear self-adjointness of perturbed 
PDEs will be used in the proof of the main results in Subsection 12.21 



Theorem 2.5 [19] If adjoint system ^ has solutions in the form v u = ef a (x,u) with 
some functions f a (x,u), then system (QJj is approximately nonlinearly self-adjoint. 

2.2 Main results 

In this subsection, we give the following discriminant criterion of approximate nonlinear 
self-adjointness for perturbed PDEs by means of nonlinear self-adjointness of the corre- 
sponding unperturbed PDEs. 

Theorem 2.6 Perturbed system (QJ) is approximately nonlinearly self-adjoint if and only 
if unperturbed system (0|) is nonlinearly self-adjoint. 

Proof. (<=). If system ([2]) is nonlinearly self-adjoint, then by Definition 12. 1} for all 
solutions u of system (j2J), there exists a substitution v a = ip a (x,u) which are solutions of 
adjoint equations (jl]), such that <p(x, u) ^ 0. Since adjoint system @ is linear in v and its 
derivatives, thus ev a = e<f a (x, u) also satisfy equations ( 
(|T]) is approximately nonlinearly self-adjoint. 



Then by Theorem I2.5[ system 

(=>). Conversely, if perturbed system ([I]) is approximately nonlinearly self-adjoint, 
then by Definition 12. 2\ there exists a nontrivial substitution w G = ^(x, u) + e<p a {x ) u) 
such that equality ( ITU]) holds. 

Observe that C = C + ew^E^, then system (JSJ) becomes 



5C 5C 6(w p El] 



5u a 5u a 



(12) 



Here, £ = w^Ep is identical to Eq.(|3]) except for the introduced variables and 
Inserting (fT2|) and w a = Wq(x,u) + ew°{x, u) into (jTUI) . that is 



5u a 



+ e 



\WQ—l[l 



5u a 



+ 



5u a 



\wi= 



+ e^)E° + e\iE}] =0(e 2 



and separating it with respect to the perturbed parameter e up to first degree, we obtain 



VU \w =i> 

6(wgEl) 8(w?E%) 

Obviously, the first equation in (|T3|) is just the required condition for nonlinear self- 
adjointness of unperturbed PDEs. 

We prove sufficiency by contradiction. Assume that system ([2]) is not nonlinearly self- 
adjoint, that means the first equation in ( !T3|) has no nontrivial solutions, i.e., ip(x,u) = 0, 
then = due to Ep ^ 0. Since 5(ip l3 Ej j )/5u ' is linear in if) and their derivatives, thus 
E\) / 8u a = and the second equation becomes 



8u a U 1= 



«■ (14) 



System (I14p has the same form as the first equation in ffTB"]) and also only has zero solutions 
4>{x,u) = 0, thus the required substitution w a = ip a (x,u) + ecj) a (x,u) = 0, that means 
system (CQ) is not approximately nonlinearly self-adjoint, which is a contradiction. Hence, 
unperturbed PDEs ([2j) must be nonlinearly self-adjoint. It proves the theorem. ■ 
Theorem E2] provides a shortcut to discriminate approximate nonlinear self-adjointness 
of perturbed PDEs. It attributes approximate nonlinear self-adjointness of perturbed 
PDEs to nonlinear self-adjointness of its corresponding unperturbed system. For the 
perturbed PDEs sharing the same unperturbed PDEs, we have the following results. 

Corollary 2.7 If the unperturbed PDEs are not nonlinearly self-adjoint, then the per- 
turbed system consisting of unperturbed PDEs with any perturbed terms are not approxi- 
mately nonlinearly self-adjoint. 

In what follows, based on Theorem 12 . 41 and Theorem [221 we gi ye a specific approximate 
conservation law formula by means of the conservation laws of unperturbed system (j2J), 
which simplify the computations of first-order approximate conservation laws of perturbed 
system (OQ). 

Theorem 2.8 Let X = X + eX\ be a first-order approximate symmetry generator of 
perturbed system (TJP, where X$ = ^d x i + n^d u a ,X\ = £{d x i + rjfd^ and the substitution 
for approximate nonlinear self-adjointness is w a = ip a + e0°"(x, u). If C = (C 1 , . . . , C n ) 
is a conserved vector of unperturbed system (TJ|) obtained by formula (TO]), then a first- 
order approximate conservation law D^iT 1 ) = 0(e 2 ) of perturbed system (TJP is given by 
the formula 

T = C\ + e + + D^W*) 5 -^ + D 3 D k (W£)^ + ... 

lv=w L ° u Suf J Su ij Su ijk 

Hic 2 + + D jW) g + WT)£| + . . (15, 



where W£ = t)q — £o M J? W[ = T)i~Ci u j an d A = w oE}i ^2 = w oEp- The Euler- Lagrange 
operators with respect to derivatives of u a are obtained from (GJ) by replacing u a by the 
corresponding derivatives, e.g. 



Proof. Observe that in conservation law formula (ITT]) , the symbols are represented as 
follows: 



thus direct computations lead to the theorem where the second-order terms of e are 



Remark 2.9 In Theorem \2.8i the conservation law Di(C l ) = is calculated via the 
substitution v a = (f a (x,u) on the solution manifold of system The approximate 

conservation law _Dj(T 4 ) = 0(e 2 ) is obtained by formula / TT3]) evaluated on the solution 
manifold of perturbed system (QP, where w a = ip a + ecj) a (x, u) and v a = <p a (x, u). 

Theorem 12.81 provides a shortcut to construct approximate conservation laws if the 
studied perturbed system consisting of two or more perturbed PDEs have the same un- 
perturbed system. It means that if unperturbed system ([2]) has a conservation law, then 
the formula (fT5l) in Theorem 12.81 gives a direct computable formula of first-order ap- 
proximate conservation laws for perturbed system (OQ). Furthermore, the approach of 
approximate symmetry definitely extends symmetry scope of perturbed PDEs [TU[T5], 
thus approximate nonlinear self-adjointness can generate new approximate conservation 
laws which cannot be obtain by nonlinear self-adjointness. 

3 Applications 

We illustrate our results by considering some examples from the literatures. Note that the 
approximate symmetries used to construct approximate conservation laws are obtained 
by the method originated from Baikov et.al [T^j . 

Example 1. Perturbed nonlinear wave equation 

The first example is to study a class of perturbed nonlinear wave equations 



where F(u) is an arbitrary smooth function. Eq. (|T6|) describes wave phenomena in shal- 
low water, long radio engineering lines and isentropic motion of a fluid in a pipe [TUE2] ■ 
In [T¥ll2"l] . two types of approximate symmetry methods were applied to Eq. (1161) and 
approximate symmetry classifications and reductions are performed. Approximate con- 
servation laws are considered via approximate nonlinear self-adjointness in [19J. 




W a = W° + eW[, JC = £ + ew?E 



omitted. 



u tt + eu t = [F(u)u x ] x , F'(u) ^ 0, 



(16) 



There exist two ways to investigate approximate nonlinear self-adjointness of Eq. (ll6p . 
The first method is based on Definition 12.21 which is performed in [19] , here, we recall 
the conclusion and use Theorem 12.61 to explore approximate nonlinear self-adjointness of 



Proposition 3.1 [7P| / Eq. I[TS\) is approximately nonlinearly self-adjoint under the sub- 
stitution 

w = (cit + c 2 )x + c 3 t + c 4 + e (^cit 2 + c 5 t + c 6 )x + ^ c 3^ 2 + c 7 t + c 8 , (17) 

where Cj, (i = 1, . . . , 8) are arbitrary constants such that w^O. 

Proof: In [6], the authors showed that the unperturbed equation of Eq. (fT6|) 

u tt - [F(u)u x ] x = 0, (18) 

is nonlinearly self-adjoint upon a substitution v = C\xt + c 2 x + c^t + C4 with arbitrary 
parameter Ci(i = 1, . . . , 4) , thus by Theorem I2.6[ Eq.( !T6|) is approximately nonlinearly 
self-adjoint. This completes the proof. ■ 
Example 2. Two perturbed wave equations having same unperturbed equation 
Consider the classical wave equation in one dimension 

u tt -u xx = 0, (19) 

which describes undamped linear waves in an isotropic medium. Obviously, Eq. f)19p is 
nonlinearly self-adjoint via a substitution v = c\u + f(x, t), where / satisfies f tt — f X x = 0. 
Hence, with the help of formula ffTTl) . we obtain conservation law formula 

C* = $>v{utt - u xx ) - v t W + vD t (W ), 

C x = Qv{u tt - u xx ) + v x W - vD x (W ), (20) 

where the infinitesimal operator X = ^d t + C,qO x + t]Qd u leaves Eq. f[T9^) invariant and, 
hereinafter, W = r) Q - £%u x - ^ u t . 

The two perturbed wave equations, which are Eq. (fT9~|) together with two perturbed 
terms, are in the form 

u u ~ u xx + eu t = 0, (21) 

and 

u u ~ u xx + e{uu t + ^tu 2 t - ^tu 2 x ) = 0, (22) 
whose approximate conservation laws are studied in [T6I22] via partial Lagrangian method. 



Since Eq. (fl9l is nonlinearly self-adjoint, thus by Theorem 12.61 both Eq.l jzTj) and (1221) are 
approximately nonlinearly self-adjoint. In particular, we have the following results. 

Proposition 3.2 1. Eq. j[2l}) is approximately nonlinearly self-adjoint upon a substitution 



w = C\U+h+e{citu+C2U+k) , where arbitrary functions h(x,t) andk(x,t) solve h tt — h xx = 
0, k tt - k xx - h t = 0. 

2. Eg. is approximately nonlinearly self-adjoint upon a substitution w = c x u + f + 



e(|citw 2 + \tfu + g) with functions f(x, t) and g(x, t) satisfying Eq. fflfy. 

Proof: We prove the second result while the first one can be done with similar proce- 
dure. 

Let the formal Lagrangian C = w[u tt — u xx + e{uu t + \tu 2 — |tw 2 )], then the adjoint 
equation is 5C/5u = Wu — w xx — e[(uw)t+t(u t w)t — t(u x w) x \. Assume that the substitution 
is in the form w = ip(x,t,u) + e<j)(x,t,u) and should satisfy equality ([TO]) , in this case, it 
becomes 



w tt - Wxx - e[(uw) t + t(u t w) t - t(u x w) x ] 

= (A + e/i )Kt - u xx ] + e\ (uu t + ^tu 2 - ^tu 2 x ), (23) 

where we omit the terms containing e 2 . 

The reckoning shows that ^(x,t,u),<p(x,t,u) and the undetermined function A ,/io 
satisfy 



VVu — ^xu — iPtu — 0, 

<f>u - A*o - *0 = 0, 2<p uu - 2ti) u - X t = 0, 

24>tu ~ i^uU + i> +ty t - X u = 0, 2<p xu - tip x = 0, (p tt 



- uiP t = 0. (24) 



Solving system (1241) gives ip = c\U+f, 4> = \c\tu 2 + \tfu+g and Ao = ci, /io = \c\tu — \tf. 

It proves the second result. ■ 
Next, we construct approximate conservation laws by means of Theorem 12.81 
Approximate conservation law for Eq. (121ft : By means of the substitution w = 

CiU+h+e(^Citu+C2U+k) and formula (fT5|) . we obtain first-order approximate conservation 

law of Eq.(EU) 



%J = C L W + e l& wu t + wW ° - w * w i + wDt{W x )] 
T M = C L W +z[t>u t + w x W x - wD x {W x ) 



(25) 



with the first-order approximate infinitesimal operator X = X + eX x = X Q + e(^\d t + 
ifd x + Vi^u) admitted by Eq.( l2~Tj) and Wi = rji — Hu x — £\u t . 

In particular, consider an approximate symmetry X = d t +e(—^ud u ), then £q = 1, £q = 
0, Wo = —u t , W\ = —hu, thus an approximate conservation law is 



T {2JJ = UtWt ~ UxxW + e 
T 6lJ = ~ u t w x + wu xt + e 



-uw t - -u t w 



—u x w uw x 

2 2 



(26) 



Induced by one approximate operator X = d t + \ud u ), formula (1261) provides more 
approximate conservations laws via substitution w = C\U + / + e(|cxtu + c 2 u + g) which 
includes the approximate conservation law in [22]. In particular, let w = \u + e^tu, then 
after proper arrangements, conservation law (1261) becomes 



1 2,1 2 

2 t 2 x 



1 

-e 

2 



tMj + uu t + tw^ 



J J2U 



1 

-WW, 

2 3 



tu x u t 



which coincides with the results in [22] up to some so-called gauge terms stated there. 

Approximate conservation law for Eq. (122ft : For Eq.f lZ?]) . with a substitution 
w = C\u + f + e(jCitu 2 + \tfu + g) with /, g satisfying Eq.f JT9|) . one has 



T) 



rpx 



Cf 



C\ 



X 

v — w 



+ e 
+ e 



CKuut + -tu 2 t - ^tu 2 x )w + (uw + tu t w)W - w t Wi + wD t (Wx] 
Cq(uih + ^tu 2 - hu 2 x )w - tu x wW + W X W! - wD x (Wij 



(27) 



with the first-order approximate infinitesimal operator X = Xq + eX\ = Xq + e(^\dt + 
i\d x + vA) admitted by Eq.f l22|) and Wi — rji — ^u x — ^\u t . 

For example, consider an approximate Lie point symmetry X = d u + e{—\tud u ), it 
means (* = = 0, W = \,W X - 



Ti 



■\tu, then by (|20|) and (1271) . we obtain 

-n't + e 



1 1 , 

j^w t u + -w(u + t«t) 

1 1 

-twu x + -tuw a 



Similarly, with substitution w = C\u + f + e(|citu 2 + \tfu + g), formula (1281) also 
gives more approximate conservations laws generated by the approximate operator X = 
d u + e(— \tud u ). Let u> = u + e^tu 2 , then fTSSl) becomes 



-«t 



-tuu t + -it' 



T J22 



u x + e 



1 



which coincides with the results in [22]. 

Both Eq. (l2Tj) and Eq. (122]) share the same unperturbed equation ffT9]) . thus formula 
( Tl5l) provides a direct and simple way to obtain approximate conservation law based on 
the known conservation laws ( 120]) of Eq. (fl9|) . Moreover, the new formulae produce new 
approximate conservation laws than the results stated in the previous literatures. 



4 Conclusion and discussion 

We propose a succinct approach to discriminate approximate nonlinear self-adjointness 
by means of nonlinear self-adjointness of the corresponding unperturbed PDEs. Conse- 
quently, a approximate conservation law formula is given based on the known conservation 
laws of the unperturbed PDEs, especially effective for the system containing two or more 
perturbed PDEs and sharing the same unperturbed PDEs. The results are applied to 



several linear and nonlinear perturbed nonlinear wave equations and approximate conser- 
vation laws are obtained. 

In addition, we consider the short pulse equation 

Uxt = u + ]:{u 3 ) xx , (29) 
6 

which was suggested in [25] as a mathematical model for the propagation of ultra-short 
light pulses in media with nonlinearities. Eq. (l29[) is not nonlinearly self-adjoint with 
substitution <p = tp(x, t,u) [BJ. Thus, by Theorem I2.6[ any perturbed equations associated 
with Eq.()29]) 

Uxt = u+ -(u 3 ) xx + ef(x,t,u), (30) 
6 

with arbitrary function f(x,t,u), is not approximately nonlinearly self-adjoint in the 
sense of Definition 12.21 However, Eq. (l29|) is nonlinearly self-adjoint under the differential 
substitution v = u t — \u 2 u x [BJ. 

As illustrated by short pulse equation, there exist some PDEs which are not approxi- 
mately nonlinearly self-adjoint in the sense of Definition 12 . 2 1 but with a general definition 
of (approximate) nonlinear self-adjointness such as differential substitution, thus it may 
consider the extended substitutions such as differential substitution containing first-order 
or higher-order derivatives. We will report these results in the forthcoming papers. 
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